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Understanding the mechanical properties of glasses remains elusive since the glass transition itself 
is not fully understood, even in well studied examples of glass formers in two dimensions. In this 
context we demonstrate here: (i) a direct evidence for a diverging length scale at the glass transition 
(ii) an identification of the glass transition with the disappearance of fluid-like regions and (iii) the 
appearance in the glass state of fluid-like regions when mechanical strain is applied. These fluid-like 
regions are associated with the onset of plasticity in the amorphous solid. The relaxation times 
which diverge upon the approach to the glass transition are related quantitatively. 
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Introduction: The understanding of the mechani- 
cal properties of amorphous solids cannot be dissociated 
from the elucidation of their structure, including the so- 
called glass transition which has eluded intensive research 
for quite some time. Plastic deformation, for example, 
arises from a qualitatively different physics in amorphous 
solids as compared to crystalline materials 0| . Since dis- 
location glide is precluded by the lack of crystalline or- 
der, much attention has been given to localized regions 
known as Shear Transformation Zones (STZ) 0, 0,3, Q 
which arc believed to be loci of plastic responses to in- 
creased strain. In this Communication we demonstrate 
that a careful identification of liquid- like structures sheds 
unexpected light on the glass transition itself, including 
a sharp definition of the transition, an identification of 
the long-sought diverging length-scale @ that accounts 
for the huge slowing down associated with the transition, 
and finally, a connection between the liquid-like defects 
and the onset of plastic deformations in the glass state. 

The system: the glass former explored here is the well 
studied system of a two-dimensional binary mixture of 
discs interacting via a softl/r 12 repulsion with a "diam- 
eter" ratio of 1.4. This two-dimensional model had been 
selected for simulation speed and, more importantly, for 
the ease of interpretation. We refer the reader to the 
extensive work done on this system 0], showing that it 
is a bona fide glass-forming liquid meeting all the crite- 
ria of a glass transition. In short, the system consists 
of an cquimolar mixture of two types of particles with 
diameter o-i = 1.4 and oi = 1, respectively, but with the 
same mass m. The three pairwise additive interactions 
are given by the purely repulsive soft-core potentials 
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where a aa = o~ a and a a b = (fa + °~b)/2- The cutoff radii 
of the interaction are set at 4.5cr a 6- The units of mass, 
length, time and temperature are m, <7i, r = o~\ \Jm/ e 
and T = e/fce, respectively, with ks being Boltzman's 
constant. A total of N = 1024 particles were enclosed in 
a square box (of area L 2 ) with periodic boundary con- 



ditions, and our simulations followed verbatim those de- 
scribed in 0], giving us a most welcome check on the 
validity of our results. While Ref. performed Molec- 
ular Dynamics, we ran both Molecular Dynamics and 
Monte Carlo simulations. The results shown below arc 
in agreement between the two simulation methods. Ref. 
0] found that for T > 0.5 the system is liquid and for 
lower temperatures dynamical relaxation slows down. A 
precise glass transition had not been identified in ■ 

Analysis: to visualize the structural features of this 
system we employ the time-honored Voronoi polygon 
construction [8(, where a polygon associated with any 
particle contains all points closest to that particle than 
to any other particle. The edges of such a polygon are the 
perpendicular bisectors of the vectors joining the central 
particle. As had been noted in the average coordi- 

nation number is 6 at all temperatures, and local coordi- 
nation numbers other than 6 are referred to as "defects" . 
Indeed, in previous work 0] the total concentration of 
"defects" was measured, sec Fig. [1] upper panel, and it 
was concluded that although there is a slight decrease 
in this concentration when slowing down sets in, noth- 
ing dramatic is observed. We propose here that a more 
careful analysis is called for, in particular to distinguish 
between "fluid-like" defects and "glass-like" defects. We 
observe that only in the liquid phase there exist small 
particles enclosed in heptagons (or even octagons), and 
large particles enclosed in pentagons (or even squares) 
(cf. Fig. [2] upper panel). In the glass phase we ob- 
serve only defects of the opposite type, i.e. small parti- 
cles in pentagons and large particles in heptagons. We 
propose that the concentration of the liquid-like defects is 
a superior indicator of the glass transition in comparison 
with relaxation times. The concentration q of these 
liquid-like defects becomes so small in the glass 
phase that we cannot distinguish it from zero (cf. 
Fig. flower panel), unless the glass is put under me- 
chanical strain, as shown below. Associated with this 
concentration we can define a typical scale, £, according 
to 
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FIG. 1: Concentration of defects under slow cooling. Upper 
panel: the average concentration of all the defects as a func- 
tion of temperature in black dots. Middle panel: the concen- 
tration of the liquid-like defects: large particles in pentagons 
(red dots) and small particles in heptagons (blue dots) . Lower 
panel: the fit of the concentration of liquid-like defects as a 
function of temperature according to Eq. @. 

Associated with the strong decrease in eg we observe a 
huge increase in the typical scale £, in agreement with 
the tremendous slowing down of the dynamics. 

In Fig. [1] middle panel wc show eg as a function of the 
temperature for a protocol of slow cooling. For tempera- 
tures larger than 0.8 the concentration follows closely an 
exponential fit, 

eg = Aexp—(AE/T) , A ks 0.094, AE ks 1.90 . (3) 

For temperatures in the range 0.3 < T < 0.8 we find an 
excellent fit to 

c.£ = B(T - T g ) 2 , 0.02, T g = 0.16 ±0.02 . (4) 




FIG. 2: Upper panel: The Voronoi polygon construction in 
the liquid state at T = 3 with the seven-color code used in 
this paper. Small particles in pentagons (heptagons) are light 
green (dark green) and large particles in pentagons (hep- 
tagons) are violet (pink). Lower panel: a typical Voronoi 
construction in the glass phase at T = 0.1. Note the total 
disappearance of liquid-like defects. 

The quality of this fit is demonstrated in Fig. [1] lower 
panel. The fit (j4|) appears to identify a sharp glass tran- 
sition T g = 0.16 ± 0.02; note, however, that there is no 
theoretical reason to expect that C£ truly vanishes at T gi 
but it becomes so small that we indeed do not see a single 
liquid- like defect in our finite-box simulations. We can- 
not exclude an exponentially small concentration that re- 
quires very much larger boxes to be observable. We can 
however state that T g is finite since we demonstrate below 
that at temperatures lower than T g the system reacts to 
shear like a solid. With the same range of temperatures 
we fit an apparently divergent length 

£(T)~ {T-T g )~ v , v = l. (5) 

The same caveat applies: in fact we can only state that 
£(T) becomes exponentially larger than the system size. 
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FIG. 3: The test of Eq. 0. The points are ln(r e ) as a 
function of £ 2 (T)/T. The solid line is a linear fit with to = 
0.71 and Afj, = 0.0055. 

At this point we connect this new structural definition 
of the glass transition to the more commonly considered 
criterion, i.e. the relaxation time (or viscosity), and show 
the one-to-one correspondence between the two. The re- 
laxation times r e for the same system were measured in 
Qvia the decorrelation times of the scattering functions, 
as a function of T, for 0.4 < T < 5 . To connect between 
the relaxation time r e and and the length scale £ we as- 
sume as usual that for the viscous fluid there exists 
a free energy of activation AG*(T) associated with the 
relaxation event, 

r e = r exp(AG*(T)/T) , (6) 

where tq is a microscopic time scale of the order of a 
single particle vibration time. The free energy of ac- 
tivation is estimated as the number of Voronoi cells 
N*(T) involved in the relaxation event, times the (tem- 
perature independent) chemical potential per cell A/i, 
AG*(T) w N*(T)Afi. The number N* is now taken as 
the typical number of Voronoi cells in regions that are free 
of liquid like defects, and therefore N*(T) w 7r£ 2 (T)/4fj, 
where f2 is the mean area of a Voronoi cell. We thus end 
up with the prediction 

r e = r exp« 2 (T)A/x/4fiT) . (7) 

Using now £ sw 1/-»/q together with the fits (jHJ) and Q 
we can compare the prediction |(7j) to the measured re- 
laxation times given in Q . This comparison is presented 
in Fig. [3] , where we see an excellent agreement with 
To = 0.71 and Afi = 0.0055. This value of Tn is of the 
order of the expected particle vibration time [31 • At the 
highest values of the temperature, T > 2, the concen- 
tration of liquid-like defects is too high, and the law (J7J) 
does not apply. 

Note that in our calculation we did not follow the tradi- 
tional steps of Adam and Gibbs who, due to the lack of in- 
dependent knowledge of N* , assume that N* /s* = N/S c 
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FIG. 4: The response of the system to shear. The measured 
energy of the system after the shear is applied is shown in 
red dots. The rhombic box is then returned to a square box 
and the energy at the end is shown in green dots. For small 
values of 8 the energy grows linearly in 5. When the value 
of S exceeds 10~ 3 the linear dependence is lost, apparently 
due to a plastic response, and upon the return to the square 
box memory is not lost and the energy does not return to the 
initial value. 



where s* is the configurational entropy of the region in- 
volved in the relaxation, whereas S c is the total config- 
urational entropy of the system. In our view this step 
is questionable, since it appears that the configurational 
entropy per particle of the whole system, including the 
liquid-like defects, is higher than the constrained config- 
urational entropy per particle of the region which is free 
of liquid-like defects. Indeed, we predict that the relax- 
ation time diverges where C£ —* 0, (i.e. at T g ) whereas 
the Adam-Gibss formula predicts divergence only when 
S c ~~ * (i.e. at the Kauzmann temperature Tfc if the lat- 
ter exists [ll| . We thus conclude that the two definitions 
of a glass transitions appear to agree on the existence of 
a finite value T g where the fluid becomes a glass. Wc 
stress that it is considerably easier to determine £(T) 
at low temperatures than to measure T e (T). The lat- 
ter quantity becomes inaccessible to measurement due to 
jamming. The concentration of liquid-like defects is much 
easier to determine since the Voronoi tcsselation remains 
highly mobile even in the glass phase. Minute changes in 
particle positions is reflected in changes in the Voronoi 
cells, allowing statistics to be accumulated at tempera- 
tures where the usual relaxation measurements become 
quite impossible. For this reason we are in a much better 
position to fit a finite value of T g . 

The onset of plasticity: having identified the glass 
transition with the disappearance of liquid-like defects, 
we proceed to strengthen this identification by analyzing 
the response of the glass state to a mechanical strain. In 
particular we show that for T <T g the system exhibits a 
finite shear modulus as is expected from a solid. To this 
aim we applied strain on the system in the glass phase, 
always starting from conditions as in Fig. [2] lower panel. 
The square box was strained to a rhombic form keeping 
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FIG. 5: A strained configuration after the linear elastic re- 
sponse is lost, 5 = 3 x 10~ 3 . Notice the reappearance of 
liquid-like defects. 



the area constant. This linear straining transformation 
is determined by the matrix h, 



h 



1 



1 + 5 



(8) 



where e is the strain tensor. The linear elastic energy 
density associated with such a transformation is 



L 2 



(2A + 4/i)(5 2 + 4n6 , (9) 



where A and \i are the 2-dimensional Lame coefficients. 
In all our simulations we applied the transformation with 
a chosen value of S, allowed the system to relax by Monte 
Carlo steps in the strained configuration, and then we ap- 
plied the inverse transformation, returning to the square 
box. The results of these simulations are most reveal- 
ing. For small values of S, 5 < ICC 3 , the system responds 
clastically, sec Fig. fH In this regime the energy increases 
linearly with <5 as expected from the leading term in Eq. 



©. As long as the response is elastic, it is also reversible: 
upon returning to the square box the energy returns back 
to its initial value. On the other hand, when the value 
of 8 exceeds 10~ 3 the linear dependence is lost, but not 
according to the elastic prediction Eq. ([9]) . Rather than 
curving up as a function of 5, the energy is curving down, 
apparently due to a plastic response; Indeed, upon the 
return to the square box memory is not lost and the 
energy does not return to the initial value. Most in- 
terestingly, for values of strain where the linear 
elastic response is lost, the strained configuration 
exhibits a significant concentration of liquid-like 
defects, see Fig. O We propose that these liquid-like 
defects are the STZ that are the loci of plastic responses. 
Quite characteristically, once these appear, the mechani- 
cal response is not reversible, as seen in Fig. 01 There is a 
remnant stress now, exhibiting the well known hystcretic 
behavior associated with the onset of plasticity. 

Summary: in summary, we have presented a novel 
view of the glass transition, identifying the transition 
as associated with the disappearance of liquid-like de- 
fects. A diverging correlation length was identified and 
measured. We demonstrated that this divergence coin- 
cides with the divergence of the relaxation time as mea- 
sured from standard time-correlation functions. The very 
same local arrangements of particles which is identified 
as liquid-like defects reappears once again when plasticity 
sets in under strain. In a subsequent publication we will 
present the statistical mechanics of this phenomenology, 
explaining and substantiating further the findings pre- 
sented here. 
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